d° Ferromagnetic Interface Between Non-magnetic Perovskites 
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We use computational and experimental methods to study d° ferromagnetism at a charge- 
imbalanced interface between two perovskites. In SrTiOs/KTaOs superlattice calculations, the 
charge imbalance introduces holes in the SrTiOs layer, inducing a d° ferromagnetic half-metallic 2D 
electron gas at the interface oxygen 2p orbitals. The charge imbalance overrides doping by vacancies 
at realistic concentrations. Varying the constituent materials shows ferromagnetism to be a gen- 
eral property of hole- type d° perovskite interfaces. Atomically sharp epitaxial d° SrTiOs /KTa03, 
SrTiOs/KNbOs and SrTiOs/NaNbOs interfaces are found to exhibit ferromagnetic hysteresis at 
room temperature. We suggest the behavior is due to high density of states and exchange coupling 
at the oxygen ti g band in comparison with the more studied d band ti g symmetry electron gas. 
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Conventionally, magnetic ordering in solids is associ- 
ated with atoms with partly occupied 3d or 5/ shells. 
There, large exchange energy J results in high-spin con- 
figuration of individual atoms. Coupling of these mag- 
netic moments depends strongly on the material in ques- 
tion. In the metallic band picture of free electrons, fer- 
romagnetic coupling is favored if the density of states at 
the Fermi level, D(Ef), is high enough. In this case, the 
exchange energy gained by spin-polarizing the carriers is 
larger than their corresponding kinetic energy increase. 
This is known as the Stoner criterion, JD(Ef) > 1. 

However, high density of states at Fermi level and high 
exchange interaction are not restricted to these shells 
only. Hund's coupling is also large in materials with 
partly filled 2p orbitals. Correspondingly, magnetic ef- 
fects in first-row elements O, N, C and B are a common 
experimental observation More specifically, holes in 
oxygen p orbitals have been widely predicted to result in 
magnetic moments , and ferromagnetic (FM) order- 
ing of these holes is obtained if the hole density is high 
enough HH-Q 

The high density of states at Fermi level is difficult to 
achieve; usually, 2p bandwidths are large, meaning the 
Stoner criterion is unlikely to be satisfied. However, in 
perovskite oxides, the HOMO oxygen ti g band is very flat 
and does not hybridize with any other bands This 
makes it an excellent candidate for realizing Stoner fer- 
romagnetism, if sufficient hole density is achieved. 

Simple charge counting dictates that combining one 
perovskite with neutral (100) layers and another with 
charged (100) layers results in an interface with addi- 
tional electrons or holes (termed n type or p type inter- 
face, respectively), to maintain charge neutrality of the 
complete structure Q. Such charge-imbalanced inter- 



faces between two perovskites have been widely studied, 
because an interface with a two-dimensional electron gas 
based on d orbitals exhibits properties such as metal- 
licity [l(j, magnetism [llj and superconductivity fl2^ . 
Magnetic ordering of these t2 g symmetry orbitals is a 
textbook example of complex behavior in the so-called 
complex oxides. 

An interface with holes on oxygen p orbitals has been 
much less studied. The oxygen 2p band will accommo- 
date holes only if both perovskites have d° occupancy in 
the bulk. One problem with the d° perovskite p type 
interface is oxygen vacancies [l3l - fl5| . which donate elec- 
trons to compensate for the charge imbalance. No 2D 
metallicity has been observed in the SrTiOa/LaAlOa p 
type interface, suggested to be either due to vacancies 
or Mott transition to an insulating state 14j | . In hole- 
doped BaTi03, however, both metallic and insulating p 
magnetic states have been obtained 

In this Letter, we show the p type interface of d° per- 
ovskites [l4, T^, 17 1 to be ferromagnetic and half- metallic. 
As the sample system, we study SrTiOa/KTaOs super- 
lattices with two similar (p or n) interfaces, starting with 
1.5 unit cells of both perovskites. This results in a struc- 
ture with one hole or one electron per supercell compared 
to nominal ionic charges, Fig. [TJ We compare the results 
of ab initio LSDA+J7, GGA+U and HSEsol hybrid func- 
tional [Hi methods, using the VASP [H] code and PAW 
[U potentials. The GGA employed is PBEsol [2l|. The 
in-plane lattice constant is fixed to the calculated SrTiOs 
lattice constant, to model epitaxy on a SrTiOs substrate. 
Ionic relaxations are performed to allow for ionic com- 
pensation of the interface dipoles, although their effect is 
expected to be very small [17| . 

Considering two-body Coulomb repulsion U or nonlo- 
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FIG. 1. (Color online) Calculated SrTiOg/KTaOg n type 
(left) and p type (right) superlattices, with nominal ionic 
charges of the layers shown. 



cal exchange (by HSEsol) is necessary due to the fail- 
ure of LSDA or GGA to properly predict the behavior 
of partly occupied oxygen p or Ti d orbitals. Magnetic 
properties are very sensitive to electron localization and 
small changes in partial occupancies, and d° magnetism 
is poorly described in local exchange-correlation approxi- 
mations [13|. For DFT+Z7, we use the Dudarev et al. ^ 
implementation. We vary the {/ ff — U — J parameter, 
the difference between local Coulomb repulsion U and 
exchange interaction J, between and 8 eV to study the 
effect of correlations. Magnetic metals yield best agree- 
ment with experiment at U values smaller than those 
calculated from first principles |24|, as the fully local- 
ized limit double counting correction exaggerates mag- 
netic moments and band gaps. Therefore, we expect U e g 
of 2 to 5 eV (on p orbitals) and 3 to 7 eV (on d orbitals) 
to yield correct behavior, based on earlier experimental 



to yield correct benavior, based on earlier experimental 
(HUE and computational @, El HE EE E3] studies con- 
sidering U and J. Further insight into suitable U e g is 
obtained by comparing the ma gne tization energies with 
HSEsol calculations, as in Ref. [28|. 

The SrTiOs/KTaOs superlattices with two p (hole- 
doped) interfaces are metallic up to U e s = 6 eV as well 
as in HSEsol. The extra holes are confined in the SrTi03 
layers only. The total p (hole-doped) superlattice mag- 
netization and energy difference between ferromagnetic 
(FM) , antiferromagnetic ( AFM) and paramagnetic (PM) 
states, as a function of U c g, are shown in Fig. [2] a). 
The DFT+t/ calculations were done in a p(2x2) super- 
cell to find possible antiferromagnetic orderings. The 
GGA/LSDA system has a partially spin-polarized FM 
ground state. A small on-site Coulomb interaction con- 
siderably increases the energy difference in favor of mag- 
netic ordering, and results in complete spin polarization 
of the conduction holes with realistic U e s values (2 to 5 
eV). With the HSEsol hybrid functional, as well, half- 
metallicity is obtained. The HSEsol FM phase energy 
is indicated in Fig. [2] a) . This would support choosing 
jTfiff = 4 eV for our system in the spirit of Hong et al. 

A charge ordered, symmetry broken AFM solution has 
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FIG. 2. (Color online) Calculated a) p type and b) n type 
SrTiOs/KTaOs 1.5-1.5 superlattice magnetization (squares, 
black) and energies of the ferromagnetic (FM, triangles, blue) 
and antiferromagnetic (AFM, upside down triangles, red) 
phases relative to paramagnetic (PM) solution, as a function 
of U e ff — U — J on oxygen 2p and Ti/Ta 3d/5d orbitals. The 
values are per p(2x2) supercell in GGA+U. The lines are 
guides for the eye. FM phase HSEsol energy and magnetiza- 
tion are marked with dashed arrows. 



2D AFM ordering of chains of parallel spins. It displays 
a Mott transition to an insulating state for U c g = 6 eV 
and higher, but it is the ground state only at unphysical 
U cS > 7 eV. 

With the n (electron-doped) superlattice, on the other 
hand, the ground state is paramagnetic metal for small 
U e g. The electrons are confined in the SrTi03 layers only. 
At large U e s and with the HSEsol hybrid functional, a 
partially spin-polarized FM state is found, but an AFM 
phase is energetically preferred (Fig. [5]b) for almost all 
U e fj. The HSEsol energy difference between the FM and 
PM phases is only 10 meV per p(2x2) supercell. This 
would suggest selecting 4 eV < U e g < 5 eV for d orbitals; 
however, as seen in Fig. [5]b), the ground state is highly 
sensitive to U c ff. The energy differences are large only 
with high U e ff, and half-metallicity or Mott insulation 
are not present. This is in accordance with paramag- 
netism observed below d 1 doping of Ti d orbitals in the 
bulk . Complex octahedral rotations fljll an d charge 
and spin ordered phases on the d orbitals [30j depend- 
ing on superlattice geometry, strain and U e g might be 
present but are not studied here, because the rotations 
are not present at room temperature in the constituent 
perovskites. 

Clearly, in the d° superlattice case, interface holes mag- 
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netize much more readily than electrons. To consider the 
effect of 2D localization on the magnetization of the elec- 
tron gas, we study a thicker superlattice, where the mag- 
netized hole or electron might spread over 5.5 unit cells 
of SrTi0 3 . The magnitude of magnetization is the same 
regardless of the thickness of the SrTi0 3 layer, demon- 
strating that ferromagnetism is a true interface effect. 
The hole density is strongly localized at the interface; in 
the middle of the 5.5 unit cell SrTi0 3 layer, i.e. only 
two unit cells away from the interface, maximum mag- 
netization density is 20 % of that at the interface. The 
strong 2D confinement is comparable to electron gas in 
n type interfaces |3lj ]. Coupling of the magnetizations 
at the two interfaces is negligible, showing that a single 
interface has a stable FM ground state. 

The presence of vacancies can be another source of 
possible magnetic signals in experiments. Ferromagnetic 
behavior in oxide materials is often attributed to oxygen 
vacancies forming localized magnetic moments, as elec- 
trons are donated to the surrounding cation d orbitals. 
In the n type interface, the lack of simple FM coupling 
in the electron-doped SrTiOs layer indicates that oxygen 
vacancies, contributing electrons to the d band, would 
not cause ferromagnetism. In the p type case, they would 
counteract hole doping and prevent magnetization. To 
study the effect vacancies have at interfaces, we calculate 
a 1.5/1.5 superlattice p(2x2) supercell with an electron 
(or hole) donating O (or K) vacancy. At intermediate U e s 
values, n type superlattices with O or K vacancies have 
close to zero magnetic moment, while p superlattices with 
O or K vacancies form completely spin-polarized metallic 
states. The calculated p superlattice magnetic moment 
is equal to the amount of holes contributed by the in- 
terface and the vacancy combined. This means that the 
interface region is dominated by intrinsic doping, and va- 
cancies will only change the total magnitude of the mo- 
ment. Metallicity is similarly obtained for 25 % oxygen 
vacancy concentration in SrTi0 3 /LaA10 3 interfaces 15| . 
indicating that a higher vacancy concentration is needed 
to change the interface to insulating. 

Finally, to find out the effect of the perovskite A and 
B ions on the localization and magnetization of the hole, 
we study similar 1.5/1.5 n and p type SrTi0 3 /KNb0 3 , 
BaTi0 3 /KNb0 3 , CaTi0 3 /KNb0 3 , SrTi0 3 /NaTa0 3 , 
BaTi0 3 /NaTa0 3 , and CaTi0 3 /NaTa0 3 superlattices, 
all of which have one non-polar and one polar dP per- 
ovskite. The obtained behavior of holes and electrons in 
all structures is identical to the model SrTi0 3 /KTa0 3 
system, which demonstrates that ferromagnetism is a 
general feature of p type interfaces of d° perovskites. 
Complete spin polarization of the holes is obtained 
with small values of £/ e ff, with energies similar to the 
SrTi0 3 /KTa0 3 case. 

To study predicted magnetization, ultrathin epitaxial 
films of d° perovskite KTa0 3 , KNb0 3 , and NaNb0 3 were 
grown onto Ti-terminated SrTi0 3 (001) single-crystal 



substrates by in situ pulsed laser deposition using a KrF 
excimer laser and high oxygen pressure (30 Pa during de- 
position and 10 5 Pa during post-deposition cooling). The 
microstructure of the KTa0 3 /SrTi0 3 interface was char- 
acterized by aberration-corrected scanning transmission 
electron microscopy (STEM) employing the high-angle 
annular dark field (HAADF) technique. Cross-sectional 
specimens for STEM analysis were prepared by stan- 
dard mechanical thinning, precision polishing and Ar ion 
milling. 
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FIG. 3. (Color online) (a) High-resolution HAADF image of 
the KTaOa/SrTi03 interface along the [100] zone axis. The 
inset shows the local fast Fourier transform pattern. The 
different elemental atoms close to the interface are marked 
by colored open circles, (b) Intensity profile of atom columns 
averaged over the indicated area from A to B in (a). 

Figure [3] a) displays an HAADF image with the view- 
ing direction along the crystallographic [100] direction of 
the SrTi0 3 substrate. From local fast Fourier transform 
(FFT) analysis (inset of Fig. [3] (a)), full epitaxial growth 
of KTa0 3 on SrTi0 3 with a KTa0 3 [100] (001) // SrTi0 3 
]100](001) relationship is derived. The growth is coher- 
ent with the in-plane lattice parameter of the KTa0 3 film 
matching that of the SrTi0 3 substrate. No dislocations 
are observed to relax the lattice strain (-2.1 %) within 
the first 10 nm of the KTa0 3 film. In high-resolution 
HAADF imaging or so-called Z-contrast imaging, the 
atoms appear with bright contrast on a dark background. 
An intensity profile across the interface (from A to B) is 
shown in Fig. [3] (b). In the KTa0 3 film, strong peaks 
from the Ta atoms are easily determined, whereas low 
intensity peaks from K atoms are not resolved. At the 
interface, two adjacent low intensity peaks due to Ti and 
K atoms are visible, indicating the formation of a p type 
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interface. Because the difference in intensity between the 
Ti and K peaks is below noise level, it is impossible to un- 
ambiguously identify the interface structure. While TiC>2 
termination of the SrTiOa substrate and the formation of 
a KO atomic layer at the heterointerface are most likely, 
some degree of intermixing and off-stoichiometry cannot 
be excluded. 

Optical constants in the KTaOs/SrTiOa heterostruc- 
tures were determined using variable angle spectroscopic 
ellipsometry and the WVASE32 software package for 
data analysis [32| ■ The model analysis considering a 
stack of substrate, film, surface roughness, and ambient 
air failed, while a nearly perfect fit was obtained when 
an additional nanolayer was introduced below the film 
(Fig. |4]a). The thickness of the KTaC>3 film determined 
from the ellipsometric data is equal to that determined 
by x-ray reflectivity, evidencing correctness of the fitting 
procedure. The thickness of the interfacial layer is about 
2 nm. The spectral features of the interfacial nanolayer 
resemble those of the SrTiC>3 substrate at large photon 
energies E > 3 eV, indicating that the origin of the in- 
terfacial layer is related to changes in electronic states 
in SrTiC>3. A profound absorption tail at E < 3 eV 
suggests the presence of in-gap states in the interfacial 
SrTiC>3 nanolayer. 




-NaNbO^ 




■KNbO, J 

- «««««3 


KTa0 3 



2 4 6 
E (eV) 



1 
H H (T) 



FIG. 4. (Color online) a) Absorption coefficient a as a func- 
tion of photon energy E in SrTi0 3 substrate (blue), KTaC>3 
film (black), and interfacial nanolayer (red). The spectra are 
obtained from ellipsometric data for the model stack shown 
schematically, b) SQUID measured total magnetization as a 
function of external magnetic field determined at room tem- 
perature in as-deposited 10x5 mm 2 films of KTa0 3 /SrTi0 3 
(red), KNb0 3 /SrTi0 3 (blue), and NaNb0 3 /SrTi0 3 (black). 



Magnetic measurements were performed in the film 
plane by a Quantum Design SQUID magnetometer 
(MPMS XL 7 T). Reciprocating sample transport en- 
abled high sensitivity (10 -8 emu). The diamagnetic re- 
sponse of the reference SrTiC>3 substrate was separately 
measured and used to extract the magnetization of the 
film and interface from the response of the film-substrate 
stack. Additionally, annealing of the samples in an oxy- 
gen atmosphere (450-500°C, 20-24 hours) was performed 
in order to check the possible influence of oxygen vacan- 
cies. Ferromagnetic-type behavior is observed in all sam- 
ples (Fig. H]b), both as-deposited and after annealing. 



The saturation moments [33| do not decrease after an- 
nealing and they remain practically unchanged with de- 
creasing temperature from 300 to 20 K, which indicates 
ordering temperatures well above room temperature. 

As stated, in d° bulk perovskites and superlattices, the 
d band is very wide, and the oxygen 2p band has a higher 
DOS and the exchange interaction J is higher for oxygen 
p than transition metal d states Q. The HSEsol hybrid 
functional yields a half-metallic density of states where 
the 2p band edges of the two spins are separated by over 
1 eV in the p interface case (Fig. [S]a)). For the n type 
interface FM phase (not pictured), on the other hand, the 
HSEsol Ti d band edge stays almost the same for both 
spins, although exchange splitting of the peak occurs. 
The strength of magnetic effects in these interfaces is 
illustrated by the HSEsol magnetization energy, which is 
340 meV for p type interface but only 10 meV for n type 
interface per p(2x2) supercell. 





FIG. 5. Densities of states for 1.5/1.5 SrTi0 3 /KTa0 3 p type 
superlattice with a) HSEsol hybrid functional, b) GGA and 
c) GGA+C/, U cS =2 eV. 

In the Stoner ferromagnetic transition, the spike in the 
DOS crosses the Fermi level when spin polarization takes 
place. This also happens in the LSDA and GGA DOS 
(Fig. [5]b)), and the system seems to be close to a fully 
half-metallic state. Half-metallic FM ground state re- 
quires a small local Coulomb repulsion (Fig. \5\ c) or in- 
clusion of nonlocal exchange (Fig. [5] a) on the oxygen 
p orbitals. The electron-doped superlattice with small 
U e g, on the other hand, retains a paramagnetic DOS vir- 
tually identical to the GGA DOS, and nonlocal exchange 
causes only partial magnetization. The hole-doped sys- 
tem clearly satisfies the Stoner criterion JD(Ef) > 1, 
while the electron-doped does not. The same mechanism 
has been suggested for hole-doped ZnO Q, where non- 
intrinsic doping methods have been considered. 

In summary, we report experimental evidence and 
ab-initio analysis of d° ferromagnetism at charge- 
imbalanced perovskite interfaces. We suggest that ferro- 
magnetic ordering and half-metallicity is present in hole- 
doped d° perovskites due to large p orbital exchange split- 
ting and the high density of states at the top of the oxy- 
gen p band under sufficient doping. The reason for the 
general scarcity of d° ferromagnets is the difficulty of dop- 
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ing the p valence band. High hole concentrations, how- 
ever, are easily obtained with intrinsic doping by inter- 
faces. Even a single interface has sufficient DOS for fer- 
romagnetism, since the carriers are strongly 2D localized 
at the interface region. Electron-doped d° perovskites 
are ordinary metals, as the bottom of the cation d band 
has larger dispersion, but they are susceptible to partial 
magnetization and possible complex magnetic ordering 
patterns under 2D localization. 
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